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Abstract 

In a previous paper we have shown how, for bosonic fields, the generating 
functional in both relativistic quantum field theory and thermal field theory 
can be evaluated by use of a standard quantum mechanical path integral. In 
this paper we extend this method to include fermionic fields. A particular 
problem is posed by Green's functions with external fermionic lines, where 
the different boundary conditions of bosons and fermions in imaginary time 
have to be accommodated within one path integral expression. The general 
procedure is worked out in the example of scalar and spinor self-energies in a 
simple model with a Yukawa coupling of a scalar to a Majorana spinor. 
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I. INTRODUCTION 



Loop corrections in thermal field theory can be calculated using a Feynman diagrammatic 
approach which is similar to the way in which analogous processes can be computed in 
ordinary quantum field theory [l]. In the Matsubara formulation of thermal field theory, 
the space is Euclidean rather than Minkowski with ( ant i-) periodic boundary conditions in 
imaginary time with extent (3 = l/(kT), where T is the temperature. Thus the frequencies 
of all fields are either integer multiples of 2tt/ (3 (for Bose particles) or half-integer multiples 
(for Fermi particles). Apart from having sums over frequencies instead of integrals, the 
perturbation theory is the same as at zero temperature. 

Recently an alternative to the Feynman diagrammatic approach has been proposed @ in 
the context of operator regularization M. In the latter, the evaluation of Green's functions 
is based on the computation of a matrix element of the form 

M xy = (x\e- im \y) (1) 

in powers of the background fields upon which the "Hamiltonian" H depends. If this is done 
by using the "Schwinger expansion" |j] 

°° roc 

e A+B = V" / daci ■ ■ ■ da n S(l — a\ — «2 . . . — a n ) 
n=i J o 

xe aiA Be a2A ---Be anA (2) 

then the calculation is quite similar in form to what is encountered in the conventional 
Feynman diagram approach; both loop-momentum and Feynman parameter integrals have 
to be evaluated. The alternative approach based on quantum mechanical path integrals 
computes Eq. (|I]) more directly, and turns out to obviate loop-momentum integrals 
completely. 

This latter method has been applied to calculations in relativistic quantum field theoret- 
ical models involving scalars, spinors, and gauge fields to one and two loop order In a 
previous paper we have shown how this approach can be used to evaluate Green's functions 
in thermal field theory in the purely bosonic case |J. Again, no loop integrals over spatial 
momenta are encountered, whereas in place of the sum over Matsubara frequencies one has 
to perform a sum over winding numbers of paths with respect to the circular imaginary 
time, the two conceptually different sums being related by Jacobi's imaginary transforma- 
tion for theta functions. The representation as a sum over winding numbers was shown to 
also constitute an interesting alternative for the derivation of high-temperature expansions 
of Green's functions. 

In this paper we consider the generalization of this alternative approach to include 
fermions, which differ in their boundary conditions in imaginary time, being antiperiodic 
rather than periodic. In the case that only fermions are involved, it is not surprising to 
find that the only difference consists of having to evaluate an alternating sum over winding 



X A similar approach has been developed independently, at one loop order, in Ref. Q. 
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numbers of paths in imaginary time, which, through the Jacobi transformation, corresponds 
to the usual sum over half-integer Matsubara frequencies. However, is not so obvious how 
the different boundary conditions have to be accounted for if both fermionic and bosonic 
fields appear in a single loop diagram, which is to be replaced by a certain quantum me- 
chanical path integral. To illustrate how this is done, we calculate the two-point functions 
at one-loop order in a simple model in which a scalar field is coupled to a Majorana spinor 
by a Yukawa interaction. The generalization to more complicated (and more interesting) 
field theories is entirely straightforward, given the methods described in Ref. [0. 



II. THERMAL GREEN'S FUNCTIONS INVOLVING FERMIONS 

We consider a four- dimensional model for a scalar field A and a Majorana spinor x with 
Lagrangian density 

1 % 1 

C = --Ap 2 A - -xi>x ~ ^X^x (3) 

(p = -id,X = X T C). 

In the imaginary-time formulation of thermal field theory the four-dimensional space 
is Euclidean, x° having the finite range < x° < (3, with Bose (Fermi) fields subject to 
periodic (anti-periodic) boundary conditions. 

We shall consider scalar and fermionic Green's functions in turn, which — in the diction 
of the background field method |7J - means that we have to introduce background fields 
for A and x-, respectively, while the other field appears only as a quantum field. 



A. Scalar Self-Energy 

With the spinor field x being a purely quantum field and A having a background con- 
tribution B, the generating functional for scalar Green's functions at one- loop order reads 

Z«[ B ]=SdeH(-f _ W ° +Afl) ) 

= const, x Det2(i|S + XB), (4) 

where "Sdet" denotes the functional superdeterminant || with the appropriate boundary 
conditions, which, for the final, ordinary functional determinant in Eq. (f|) are antiperiodic 
ones, because the latter is the fermionic part of the former. 
Introducing the matrix operator 

+ XB) 2 = p 2 - iX^B + Bj>) - X 2 B 2 
= (p - iX-iBf - 3X 2 B 2 , (5) 

we regulate Eq. |] by the £ function [^,0 

C(s) = Tr— - dtt s -\' m (6) 
I s Jo 
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according to In [B] = jC(0). We therefore have to evaluate 

Tr exp(— Ht) = J d 4 xd 4 y 8{x — y) (x\tr exp(— Ht)\y). (7) 

While at one-loop order only the diagonal matrix elements enter, at higher loop orders the 
general matrix elements are the basic building blocks of the generating functional, if defined 
by the method of operator regularization [|J. 

In order to obtain the one-loop contribution to the (thermal) scalar self-energy, only the 
term bilinear in the background field B needs to be extracted from Eq. (]?]). 



1. Evaluation by Schwinger expansion 

We first evaluate Eq. (|7|) by employing the Schwinger expansion, Eq. (g), which re- 
produces the expressions derived from the more customary Feynman rules together with 
Feynman parametrization of denominators. This will serve as a reference for the results to 
be derived from the use of quantum mechanical path integrals. 

To second order in B, Eq. (Q) yields 



Tr exp(-Ht) 



A Tr 



t B 2 e- p2t - - I 1 da e- {l - a)p2 \i)B + B^)e~ apH ^B + Bj>) 
2 jo 



(8) 



Evaluating further by insertion momentum states \p), we have, remembering the an- 
tiperiodic boundary conditions for the functional determinant in Eq. (f|) and hence for the 
functional trace in the subsequent expressions, 



(x\p) 



1 



(27T) 3 / 2 v^ 



exp % 



2ir 



1 , 



and 



1 2tt 
{plBlq) = (27r) 3 / 2 ^ jB(P ~ q ' J {Hp ~ Uq))) 



(9a) 



(9b) 



with n p and n q integers. After appropriate shifts of the integration and summation variables, 
and then using / e~ q t d 3 q/(2n) 3 = (47rf)~ 3//2 , we end up with 



Tr exp(-Ht) 



. 2tt , 
xB(p,—n p ) 



-4A 2 



B 2 



d 3 p 



■i 



dae~ a{1 ~ a)p l a ( n p+ n i+V 2 +(i-")K+|) 2 ]t 

2tt 



^ - a(l - 2a)p 2 + (j) 2 (n p (n q + + 2(n q + ^) 2 



B(-p, 



n p ). (10) 



In accordance with the Feynman rules in the Matsubara formalism, we encounter imaginary 
frequencies for the scalar background fields which are integer multiples of 2tt//3, and there 
is a further summation over half-integer (n q + |) multiples corresponding to the frequencies 
of the internal fermion lines. 
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2. Evaluation by a quantum mechanical path integral 



We now compute the matrix element appearing in Eq. (|?p by using the representation 
through the standard quantum mechanical path integral pll 



( x | exp |_[I(p_^)2 + V ] t | 

IP J Dqexp J dr 



~q 2 (r)+iq(r)-A(q(r))-V(q(r) 



(11) 



with g(0) = y and q(t) = x. IP stands for path ordering in the case that A and V are matrix 
valued. 

Comparing with Eq. fl5|) and rescaling t — > t/2 for convenience, we identify A = iX'jB 



and V = |A 2 -B 2 . Consequently we have 



Tr exp(-Ht) = TrW J Dqexp J dr 
To second order in B this becomes 



Tr exp(-Ht) 



- l - q \T)- q {T). lB -\\*B> 



(12) 



B 2 



A 2 Tr J Dqe^Io dT ^(-^J o dr.B 2 ^)) 

+ f dr x r dr 2 g(n) ■ jB(q(n)) q{r 2 ) ■ ~fB{q{r 2 )) 



(13) 



where in the last term the integration over r 2 is restricted by the necessity of path ordering 
of the 7 matrices.] 2 ] 

If we now take B to be a plane wave field, 



B{q{n)) 



(27r) 3 / 2 v^ 



(27T) 3 / 2 v^ 



3 j[k l q(r i )+(^)n i( j4(Ti)] 



(14) 



and assign labels 1 and 2 to the 7 matrices contracted with q{r{) and q(r 2 ), respectively, 
then we can write 



Tr exp(-Ht) 



\%b 2 
~ (2vr) 3 /3 



Tr j Dqe-*&**VL~ f dr x e/o *-fo(rM*i+W 



_|_ f ^ ^ T ^ e J^dT[q(T)iik 1 8(T-T 1 ) + ik 2 S(T~T 2 )-'-/lS'(T-T 1 )-^ 2 8'{T-T 2 ))} 



provided we keep only the linear terms in 71 and 72. 
The path integral is now of the form 



2 An alternative to explicit path ordering of 7 matrices, which moreover removes their explicit 
appearance altogether, could be given by the introduction of fermionic path integration variables 
in addition to q [11|. However, this has no bearing on the problem of boundary conditions in 
imaginary time, which is our main concern here. 
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-_g(r)+ g (r).r(r) 



(16) 



where the superscript /3 reminds of the fact that g 4 lives on a circle with circumference f3, i.e. 
g 4 = g 4 + n/3, with n integer. This topological constraint prevents us from a direct evaluation 
by a shift of variable 



5(r) = g(r)- f dr'G{r^')T{r') 
Jo 



as in zero-temperature field theory [0. In the latter, 



- t — r 

2 1 



- r + r 
2 V 



rr 



(17) 



(18) 



would render the path integral purely Gaussian. 

The problem is that of a path integral quantization of a particle on a circle, and the 
general solution has been given in Ref. \T1\ in terms of a superposition of unconstrained 
path integrals, 



(19) 



which corresponds to summing over all paths from y to x which have different winding 
number in the fourth variable. 

In the purely bosonic case treated in Ref. ||, the phases 5 n could be set to zero. The 
present path integral, however, is supposed to give the matrix elements for taking a functional 
trace with antiperiodic boundary conditions in imaginary time. We therefore have to require 
that P(L. y ) picks up a phase equal to —1 when making one turn in the circular imaginary time. 
This determines S n = nir, and P@ is given by an alternating sum over winding numbers. 

On the left hand side of Eq. (|l^), we can indeed perform the shift of Eq. (|17|), yielding 



poo 

fay) 



[r] 



(2^F exp 



(x - yf 



2t 



t Jo 



dr[xr + y(t — r)] • T(r) 



~ drdr'G(r, r')T(r) ■ r(r') 
2 Jo 



(20) 



with G given by Eq. p8]) . 

Evaluating the above path integral in this manner and computing the trace by integrating 
over x after setting x = y, which has the effect of equating k\ = —k 2 = k, we obtain 



Tr exp(-Ht) 



B 2 



2A 2 M 2 
P(2xt) 2 



^(-l)» e -«W*)|t 

rt *• 



+ f dn dT 2 e -k 2 *(l-")t-in/3k4a 

Jo Jo 



i-k>-i) 2 + (f -^ 4 (a-i)) 2 



(21) 



where a = {r\ — T 2 )/t. Because fc 4 = jj-nk with integer nk, we find that the integrand f(a) 
in Eq. ( PT|) is symmetric under a <-> 1 — cc together with n <-» — n, which allows us to simplify 
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f dn r dr 2 f(a) - t 2 f 1 da(l - a)f(a) = % C daf(a). (22) 
Jo J o jo 2 Jo 

Eq. pi|) involves an infinite sum different from the one encountered in the result of the 
more conventional Schwinger expansion, Eq. (|l0|). Introducing Jacobi's 9 functions 

oo y 

6 2 (z\r)= J2 exp(tTT[r(n + -) 2 + 2z(n + -)]), (23a) 

n=— oo 

oo 

exp (ivr[™ 2 + 2zn\) , (23b) 



n=— oo 
oo 



4 (z|r) = ^ (-l) n exp (iir[Tn 2 + 2zn\) , (23c) 



(25) 



the path integral result Eq. (BTJ) can be expressed in terms of 84 and its derivatives with 
respect to its two arguments, whereas the conventional result Eq. fllCf ) requires Q%. Jacobi's 
so-called imaginary transformation 

9 3+i (z\r) = (- lT )-\e-^ 2 ^9,J- - -), (24) 

where i = —1,0, or 1, relates #4 to 9 2 , and indeed, after some algebra, establishes the 
equivalence of the two results. 

B. Spinor Self-Energy 

We shall now consider the case of Green's functions with external fermionic lines, which 
receive radiative corrections from both (quantum) scalar and fermionic degrees of freedom 
at the same time. By taking now A to be a purely quantum field and \ having a background 
part 7], the generating functional for fermionic Green's functions at one-loop order reads 

p 2 Xf] \ 
Xrj ip 1 J ' 

Dropping an irrelevant constant over-all factor, we may rescale the supermatrix Y in 
Eq. ( p5|) with the field-independent supermatrix 

*-(!-;)■ (26) 

and regulate the superdeterminant of H = XY using the ( function 

I r°° 

C(s) = Str— - dtt s ~ l e- Ht (27) 
T(s) Jo 

according to lnZ^ 1 ^?]] = — 1^'(0). We thus are confronted with having to evaluate 

where "Str" denotes the supertrace including a functional trace with the appropriate bound- 
ary conditions and "str" the supertrace in the purely algebraic sense. 

In order to obtain the one-loop contribution to the (thermal) spinor self-energy, we now 
have to extract the term bilinear in rj from Eq. (p8|) . 



Z (1) [q\ = Sdet-5Y[7?] = SdeH ( 
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1. Evaluation by Schwinger expansion 



Again, we first lay out the derivation of the thermal spinor self-energy employing the 
Schwinger expansion, Eq. (0), by which we formally have 



(x|strexp(— Ht)|y) 



</'/ 



i\H 2 f 1 c/«(x|e- (1 - a)p2 '^e- Qp2 'r/|y). 
Jo 



(29) 



In order to correctly take into account the (anti-)periodic boundary conditions of the su- 
permatrix operator H when computing the functional trace in Eq. (28) in momentum-space, 
one has to remember that fj appears in its Bose-Fermi block, and 77 in the Fermi-Bose one. 
Therefore, 



(pM?) 



1 , 2tt 2tt. 1 



(30) 



with n p and n„ integers, which leaves us with 



Strexp(-Ht) 



-i\H 2 I da 



nn 



d 3 p d 3 q 



J (2tt) 3 /^ 



E 



-[(l- a )(p 2 +(27rn p //3) 2 +a(q 2 +(27r(n 9 + i)//3) 2 )]t 



xmp-q,—{n p -n q - -)) 



2tt 



1 27T . 1 , . 

^(q-p,-^K + - -np)). (3i) 



With the shifts p — > p + q, n p — > n p + n q , followed by q — > q — ap, and then using the 
integral / e - q2 *c/ 3 g/(2vr) 3 = (4vrt)- 3 / 2 , Eq. (0) reduces to 



Strexp(-Ht) 



iA 2 t 2 / cfo 



'/'/ 



d 3 p 



(4^)3/2/3^ 



-{a(l-a)p 2 +(2ir/0) 2 {(l-a)(n p +n q )2+a(n q + ±) 2 }]t 



2vr 



xr/(p, — K--)) 



2vr 



-(l-a)p / + — K + -) T 



, 27T . 1.. 

77(-p,--g-(n p --)). (32) 



Evaluation by a quantum mechanical path integral 



In order to represent Eq. (|28|) by the standard quantum mechanical path integral Eq. (|TT]) . 
we rewrite 



H 



p 2 Xrj 
-iXp'r/ p 2 



Sr + 

¥7^ /V / 



Xfj\ 



(33) 



leading to 



(x\exp(-Ht/2)\y) 



P / Dqexp J dr 



1 . 2 .. / 
2 9 + ^U^0 







-77^ 



(34) 
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where we have again rescaled t — > t/2 in Eq. fl28|) in order to employ the standard quantum 



mechanical result Eq. (11) directly. 

To second order in rj, fj we find from Eq. (|34| ) that 



(a?|strexp(-Ht/2)|j/) 
J Dqe^fo^^sti 



vv 



2! 



p n / ^ 



-Mm + huVi,u o 



(35) 



where q { = q(r^ and r] { = r)(q(r^). 

Computing the (algebraic) supertrace gives two contributions, the trace over the bosonic 
part minus the one over the fermionic part, which are almost identical apart from the sign, 
were it not for the path ordering prescription IP, which is required because we are dealing 
with a matrix valued "Hamiltonian" H. The bosonic part of the supertrace reads 



(x|tr B exp(-Ht)|y) 

vv 



2»7l 



hV2 ~ lpuV2,» , ■ 



(36) 



whereas the fermionic part equals minus this expression where the indices 1 and 2 are 
interchanged except in the integration symbols. 

Without the topological constraint of a circular imaginary time q^, the bosonic and the 
fermionic trace could be simply combined, resulting in a symmetrized integration domain 
in n and Ti. However, in the preceding section we have seen that the path integral corre- 
sponding to a fermionic trace when expressed in terms of an unconstrained path integral, 
acquires a nontrivial phase for odd winding numbers in imaginary time, whereas the bosonic 
part does not. The two pieces have therefore to be treated as two different path integral 
expressions in the thermal case. 

Keeping this difference in mind, we can transform Eq. ([36|) and its fermionic counterpart 
into a Gaussian path integral by taking the background wave functions to be plane wave 
fields whose frequencies are half integral multiples of 2ti//3, 



This reduces Eq. 



Ui 



(2^)3/2^ 



(2^)3/2^ 



3 4 k i<lfo)+(^r)(ni+§)g4(ri) 



(37) 



to 



(x|tr B exp(-Ht)|2/) 
A 2 



'/'/ 



4(2tt) 3 /3 Jo Jo 



1 dri [ T1 dr 2 u 1 [ Dqfefo ^[-59 2 ( T )+9(' r )-(-* fe i <5 ( r - T i)+^25(r-^)-7(5'(T-n 2 ))] 



.LU e Io dT[-\q 2 {T)+q{T)-(-ik 1 5(T-T 1 )+ik 2 S{T-T2))\ , ? ! 

2 J lin. 



(38) 



where the subscript "lin." indicates that we have exponentiated 7 matrices with the provision 
that only first order terms in 7 are to be kept. 
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Eq. (pq) is now of the form of Eq. (|T6| ) with linear sources T = —iki6(T — Ti)+ik2S(r — T2) 
in the first term, and T' = T — j8'(t — r 2 ) in the second one. The fermionic part of the 
supertrace gives rise to a similar expression with source terms in which the labels 1 and 2 
are interchanged. The two path integral contributions can now be evaluated by Eq. ( |20|) 
and Eq. (0), taking into account the alternating signs picked up by the one corresponding 
to the fermionic part of the supertrace. Performing also the integration over x after setting 
x = y, which enforces momentum conservation so that k\ = k2 = k, we arrive at the result 



A 2 



Strexp(-Ht/2) = 

u(-k){e- inf3k4a 



1 fcri n 2 (j 2 , 



3 n/3 74 /t+f/(2a-l) _ ^ 



lin. 



u(k), 



(39) 



with a = <7\ — o"2 and Oi = Ti/t. 

Because k^ = y-(n,k + |), we find that the integrand of Eq. (^) is symmetric under 
a <-> 1 — a, since then £- %n P k ± a (-l) n e m,3<;4a . Therefore we can replace Jq 1 do\ J CT1 <icr2 - ^ 
| Jq da. The same substitution, when applied to only part of the integrand, allows us to 
write Eq. (|39D either as a straight or an alternating sum over winding numbers n, which can 
be written in terms of Jacobi 9 functions, Eq. (|2"B|), as 



Strexp(-Ht/2) 
A 2 n 



4(2?rt)- Jo 



4(2vrt) 2 Jo 



dae~ k a{1 - a)t u(-k) 
dae-^ k2a(1 - a)t u(-k) 



03 (" 

04 (" 



+ ^7i ^^ p -»/(l- Q ) 



2tt 



2nit 



2ixt> 



u{k) 



fik^a /3j4 i/3 



2tt 



2nit 



- lin. 



(40a) 
(40b) 



lin. 



where because only the terms linear in 7 matrices are to be kept we can treat them as if they 
were c-numbers. Each of these comparatively compact expressions can easily be shown to be 
equivalent to the standard result Eq. fl52|) by means of the Jacobi transformation, Eq. ( p4| ) . 
The two representations of Eq. (|40|) are thus related to sums over integral and half- integral 
Matsubara frequencies, respectively, corresponding to the two possible ways of reading the 
standard expression primarily as a sum over the internal bosonic or fermionic frequencies. 



III. SUMMARY 

We have shown that the alternative method of evaluation of Green's functions based on 
the standard formulae for quantum mechanical path integrals can be adapted to the case of 
thermal Green's functions involving fermions, provided that a separation is made between 
purely bosonic and purely fermionic contributions. This separation is a natural one when 
the generating functional is formulated in terms of supermatrices. The resulting quantum 
mechanical path integrals are topologically constrained, the one associated with fermionic 
contributions picking up a minus sign for increasing the winding number by one. The 
evaluation of Green's functions along these lines leads to Feynman-parametrized results, the 
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Feynman parameters being related to the points of insertion of external lines, but no loop- 
momentum integral is encountered. In place of the usual sum over Matsubara frequencies, 
however, a sum over winding numbers emerges. The latter has been shown in Ref. || to be 
more directly amenable to a high-temperature expansion. The generalization of the results 
presented here for a simple four- dimensional model of scalars with Yukawa coupling to 
fermions to more complicated ones can performed with the techniques developed in Ref. @ . 
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